Initial eccentricity fluctuations and their relation to higher-order flow harmonics 
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Monte Carlo simulations are used to compute the centrality dependence of the participant eccen- 
tricities (s n ) in Au+Au collisions, for the two primary models currently employed for eccentricity 
estimates - the Glauber and the factorized Kharzeev-Levin-Nardi (fKLN) models. They suggest 
specific testable predictions for the magnitude and centrality dependence of the flow coefficients v n , 
respectively measured relative to the event planes They also indicate that the ratios of several 
of these coefficients may provide an additional constraint for distinguishing between the models. 
Such a constraint could be important for a more precise determination of the specific viscosity of 
the matter produced in heavy ion collisions. 
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Collective flow continues to play a central role in on- 
going efforts to characterize the transport properties of 
the strongly interacting matter produced in heavy ion 
collisions at the Relativistic Heavy Ion Collider (RHIC) 
(lHTfjj . An experimental manifestation of this flow is the 
anisotropic emission of particles in the plane transverse 
to the beam direction [I?], HH. This anisotropy can be 
characterized by the even order Fourier coefficients; 
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where <f> p is the azimuthal angle of an emitted particle, 
^>rp is the azimuth of the reaction plane and the brackets 
denote averaging over particles and events . Charac- 
terization has also been made via the pair-wise distribu- 
tion in the azimuthal angle difference (A<fi = (f>i — fa) 
between particles [H IM 1J ; 
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Anisotropic flow is understood to result from an asym- 
metric hydrodynamic-likc expansion of the medium pro- 
duced by the two colliding nuclei. That is, the spacial 
asymmetry of the produced medium drives uneven pres- 
sure gradients in- and out of the reaction plane and hence, 
a momentum anisotropy of the particles emitted about 
this plane. This mechanistic picture is well supported by 
the observation that the measured anisotropy for hadron 
Pt < 2 GeV/c, can be described by relativistic hydrody- 
namics @, M M, 0, EE HM2 • 

The differential Fourier coefficients v 2 (N part ) and 
V2{pt) have been extensively studied in Au+Au collisions 



22 at RHIC [2CJ, l32H38j . One reason for this has been the 

23 realization that these elliptic flow coefficients are sensi- 

24 tivc to various transport properties of the expanding hot 

2 5 medium 0-Q, i, GH, [H, [si . Indeed, considerable 

26 effort has been, and is being devoted to the quantita- 

27 tive extraction of the specific shear viscosity 77/ s (i.e. the 

28 ratio of shear viscosity 77 to entropy density s) via com- 

29 parisons to viscous relativistic hydrodynamic simulations 

30 [9l-[T2l. [Til. [T5I |30| . transport model calculations [rl [l3l. |42| 

31 and hybrid approaches which involve the paramctriza- 

32 tion of scaling violations to ideal hydrodynamic behavior 

33 0(11,11(1 El, El- The initial eccentricity of the collision 

34 zone and its associated fluctuations, has proven to be an 

35 essential ingredient for these extractions. 

36 Experimental measurements of the eccentricity have 

37 not been possible to date. Consequently, much reliance 

38 has been placed on the theoretical estimates obtained 

39 from the overlap geometry of the collision zone, speci- 

40 fied by the impact parameter b or the number of par- 

41 ticipants iVp ar t [3ll. l34l. l43l. I45T452T ] . For these estimates, 

42 the geometric fluctuations associated with the positions 

43 of the nuclcons in the collision zone, serve as the under- 

44 lying cause of the initial eccentricity fluctuations. That 

45 is, the fluctuations of the positions of the nucleons lead 

46 to fluctuations of the so-called participant plane (from 

47 one event to another) which result in larger values for 

48 the eccentricities (e) referenced to this plane. 

49 The magnitude of these fluctuations arc of course 

50 model dependent, and this leads to different predictions 

51 for the magnitude of the eccentricity. More specifically, 

52 the £2 values obtained from the Glauber 

53 factorized Kharzeev-Levin-Nardi (fKLN) 




531 and the 
55j models, 
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(the two primary models currently employed for eccen- 
tricity estimates) give results which differ by as much 
as ~ 25% [H, [53] - a difference which leads to an ap- 
proximate factor of two uncertainty in the extracted rj/s 
value [1, HI]. Thus, a more precise extraction of rj/s re- 



2 



: quires a clever experimental technique which can measure 
i the eccentricity and/or the development of experimental 
i constraints which can facilitate the requisite distinction 
; between the models used to calculate eccentricity. 

Recently, significant attention has been given to the 

■ study of the influence of initial geometry fluctuations on 
i higher order eccentricities £„ in >3 HH, E3, I5014521 [58l — 
1 [60[ , with an eye toward a better understanding of how 
i such fluctuations manifest into the harmonic flow corre- 
lations characterized by v n (for odd and even n), and 

i whether they can yield constraints that could serve to 
: pin down the "correct" model for eccentricity determina- 
i tion. For the latter, the magnitude of e n and its detailed 
i centrality dependence is critical. Therefore, it is essen- 
i tial to resolve the substantial differences in the e n values 

■ reported and used by different authors [3(], HH, [50l — 



Here, we argue that the magnitudes and trends for 
the eccentricities e„ imply specific testable predictions 
for the magnitude and centrality dependence of the flow 
coefficients v n , measured relative to their respective event 
planes *& n . We also show that the values for e n obtained 
for the Glauber [H, H and fKLN gl |H| models, in- 
dicate sizable model dependent differences which could 
manifest into experimentally detectable differences in the 
centrality dependence of the ratios W3/(«2) 3 ^ 2 , va/{v2) 2 
and V2/v n ^ n >3- Such a constraint could be important for 
a more precise determination of the specific viscosity of 
the hot and dense matter produced in heavy ion colli- 



I. ECCENTRICITY SIMULATIONS 

Monte Carlo (MC) simulations were used to calcu- 
late event averaged eccentricities (denoted here as e n ) in 
Au+Au collisions, within the framework of the Glauber 
(MC-Glauber) and fKLN (MC-KLN) models. For each 
event, the spatial distribution of nuclcons in the collid- 
ing nuclei were generated according to the Woods-Saxon 
function: 



p(r) = 



Po 



1 + e (r-flo)/d' 



(3) 



where Rq = 6.38 fm is the radius of the Au nucleus and 
d = 0.53 fm is the diffuscness parameter. 

For each collision, the values for iV part and the num- 
ber of binary collisions N co n were determined within the 
Glauber ansatz (53[. The associated e n values were then 
evaluated from the two-dimensional profile of the density 
of sources in the transverse plane p s (r±), using modified 
versions of MC-Glauber M and MC-KLN [55j respec- 
tively. 

For each event, we compute an event shape vector S n 
and the azimuth of the the rotation angle for n-th 



harmonic of the shape profile [13, HO] ; 

S„x = 5„cos(n*„) = J dr_ L/ o s (r_ L )w(rj_)cos(n<^t) 
S ny = S n sm{n^ n ) = J dr±p 8 (rj_)u)(r_ L )sin(n 



tan 



(6) 



42 where (j) is the azimuthal angle of each source and the 

43 weight w(rj_) = r± 2 and w(r_i_) = r± n are used in 

44 respective calculations. Here, it is important to note 

45 that the substantial differences reported for e n in Rcfs. 

4 6 [U El, S3, \EM$ is largely due to the value of 

47 oj(t±) employed. 

The eccentricities were calculated as: 



and 



e n = (cosn(0 - *„)) 



(cosn(0 — 'I'm)) , nj^m. 



(7) 
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48 where the brackets denote averaging over sources and 

49 events belonging to a particular centrality or impact pa- 

50 rameter range; the starred notation is used here to dis- 

51 tinguish the n-th order moments obtained relative to an 

52 event plane of a different order 

For the MC-Glauber calculations, an additional en- 
tropy density weight was applied reflecting the combina- 
tion of spatial coordinates of participating nucleons and 
binary collisions [H, [56[ ; 
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oc 
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53 where a 

54 ments as a function of iV par t for Au 



0.14 was constrained by multiplicity measure- 

Au collisions [oil ]. 

55 These procedures take account of the eccentricity fluctu- 

56 ations which stem from the event- by- event misalignment 

57 between the short axis of the "almond-shaped" collision 

58 zone and the impact parameter. Note that e n (cf. Eq. [7]) 

59 corresponds to v n measurements relative to the so-called 
eo participant planes [34L [53l |. That is, each harmonic e n 
ei is evaluated relative to the principal axis determined by 

62 maximizing the n-th moment. This is analogous to the 

63 measurement of v n with respect to the n-th order event- 

64 plane in actual experiments [62]. It however, contrasts 

65 recent experimental measurements in which a higher or- 

66 der coefficient (V4) has been measured with respect to a 
e? lower order event plane (^2) [H,[6^]. Note as well that 

68 we have established that the angles ^ n for the odd and 

69 even harmonics are essentially uncorrelated for the -/V par t 
7» range of interest to this study. 



A. Results for ui(r±) = r± 2 and w(rj_) = ri" 

Figure [1] shows a comparison of £„.„<6 vs - IVpart f° r 
u(r±) = r ± 2 , for MC-Glauber (a) and MC-KLN (b) for 
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FIG. 1. Calculated values of e„,„<6 vs. AT par t for ui(r±) = ri 2 
for MC-Glauber (a) and MC-KLN (b) for Au+Au collisions. 
The open and filled symbols indicate the results for odd and 
even harmonics respectively. 



75 Au+Au collisions. The filled and open symbols indi- 

1 cate the results for the even and odd harmonics respec- 

2 tively. For this weighting scheme, e„ is essentially the 

3 same for n > 3, and have magnitudes which are signifi- 

4 cantly less than that for £2, except in very central colli- 

5 sions where the effects of fluctuation dominate the mag- 
e nitude of e n . n >2- Note the approximate l/\/(-/Vp ar t) de- 
7 pendence for e n ,n>3- The smaller magnitudes for e n , n >3 
a (with larger spread) apparent in Fig. D^b), can be at- 

9 tributed to the sharper transverse density distributions 

10 for MC-KLN. 

11 Figure[2]shows a similar comparison of e n ,n<6 vs. N par t 

12 for calculations performed with the weight cj(t±) = r± n . 

13 This weighting results in an increase in the sensitivity 

14 to the outer regions of the transverse density distribu- 

15 tions. Consequently, the overall magnitudes for e Iljrl >3 
is are larger than those shown in Fig. [TJ This weighting 
17 also lead to a striking difference in the relative magni- 
la tudes of £„,„> 2 for MC-Glauber (a), MC-KLN (b) and 
19 the results for ui(r±) = r± 2 shown in Fig. [TJ 




FIG. 2. Same as Fig. [T]for w(r x ) = r ± n . 

29 catc that rj/s is small 0, @, 0, IMl, M, \M El H - of 

30 the same magnitude as for the conjectured KSS bound 

31 77/s = 1/(4tt) [H. 
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20 II. ECCENTRICITY RATIOS 

21 The magnitudes and trends of the calculated eccentric- 

22 ities shown in Figs. [1] and [5] are expected to influence the 

23 measured values of v n . To estimate this influence, we 

24 first assume that the resulting anisotropic flow is directly 

25 proportional to the initial eccentricity, as predicted by 

26 perfect fluid hydrodynamics. Here, our tacit assumption 

27 is that a possible influence from the effects of a finite 

28 viscosity (r]/s) is small because current estimates indi- 



FIG. 3. Comparison of £2,4,5/£3 vs. Apart for Au+Au colli- 
sions. Results are shown for MC-Glauber (a) and MC-KLN 
(b) calculations. 

Figure [TJ indicates specific testable predictions for the 
relative influence of e n , n >2 on the magnitudes of v n , n >2- 
That is, (i) £2 should have a greater influence than £„ irl >3 
in non-central collisions, (ii) the respective influence of 
£n,n>3 on the values for w„. n >3 should be similar irre- 
spective of centrality and (iii) the ratios f 4, 5,6/^3 should 
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FIG. 4. Same as Fig. ED for u>(r±) = r ± " . 

38 follow a specific centrality dependence due to the influ- 

1 ence of £4,5,6/£3- Such a dependence is illustrated in Fig. 

2 E2 where we show the centrality dependence of the ratios 

3 £2,4, 5/ S3, obtained for MC-Glauber (a) and MC-KLN (b) 

4 calculations. They suggest that, if MC-Glauber-like ec- 

5 centricities, with weight u>(r±) = r± 2 , are the relevant 
e eccentricities for Au+Au collisions, then the measured 
7 ratio V2/V3 should increase by a factor m 2, from cen- 
s tral to mid-central collisions (A part ~ 350 — 150). For 
Ap ar t % 150, Fig. EJa) shows that the ratio V2/V3 could 

10 even show a modest decrease. The eccentricity ratios in- 

11 volving the higher harmonics suggest that, if they are 

12 valid, the measured values of ^4,5,6/^3 should show lit- 
is tie, if any, dependence on centrality, irrespective of their 

14 magnitudes. 

15 The ratios £2,4,5/^3 obtained for MC-KLN calculations 

16 are shown in Fig. EH (b). While they indicate qualita- 

17 tive trends which are similar to the ones observed in Fig. 
is EH (a), their magnitudes and their detailed dependence 

19 on centrality are different. Therefore, if the qualitative 

20 trends discussed earlier were indeed found in data, then 

21 these differences suggest that precision measurements of 

22 the centrality dependence of the relative ratios for V2/V3, 

23 V4/V3, 115/113, ... for several pr selections, could provide a 

24 constraint for aiding the distinction between fKLN-likc 

25 and Glauber-like initial collision geometries. Specifically, 

26 smaller (larger) values of the relative ratios are to be ex- 

27 pected for V2/V3 and V4/V3 for Glauber-like (fKLN-like) 

28 initial geometries. Note the differences in the expected 

39 centrality dependencies as well. 

31 Figure 2| compares the eccentricity ratios £2, 4, 5/^3 ob- 

32 tained for MC-Glauber (a) and MC-KLN (b) calcula- 

33 tions with the weight w(r±) = r± n . The magnitudes 

34 of these ratios and their centrality dependencies are dis- 

35 tinct for MC-Glaber and MC-KLN. They are also quite 



FIG. 5. Comparison of Ez/{E2) A ^ 2 vs. iVpart (a) and ea/{£2) 2 
vs. iVpart (b) for MC-Glauber and MC-KLN initial geometries 
(as indicated) for Au+Au collisions. 



36 different from the ratios shown in Fig. E3 This suggests 

37 that precision measurements of the centrality dependence 
3s of the relative ratios V2/V3, V4IV3, V5/V3,... (for several 

39 pt selections) should not only allow a clear distinction 

40 between MC-Glauber and MC-KLN initial geometries, 

41 but also a distinction between the the w(rj_) = r± 2 and 

42 w(rj_) = r_i_" weighting methods. 

43 A finite viscosity will influence the magnitudes of v n . 

44 Thus, for a given p T selection, the measured ratios for 

45 V2/V3, V4/V3, V5/V3, ... will be different from the eccen- 

46 tricity ratios shown in Figs. [3] and [4] Note as well that, 

47 even for ideal hydrodynamics, the predicted magnitude 

48 of «4/£4 is only a half of that for V2/S2 [1^. Nonethe- 

49 less, the rather distinct centrality dependent eccentricity 

50 patterns exhibited in Figs. E3 and [4] suggests that mea- 

51 surements of the ratios of these flow harmonics should 

52 still allow a distinction between MC-Glauber and MC- 

53 KLN initial geometries, as well as a distinction between 
5s the two weighting methods. 

56 The ratios ^/(t^) 3 / 2 and Vi/fa) 2 have been recently 

57 found to scale with pt [H[ , suggesting a reduction in the 

58 influence of viscosity on them. Thus, the measured ratios 
» v n /{v 2 ) n/2 could give a more direct indication of the cen- 

60 trality dependent influence of £ n /(£2)"^ 2 on v n / (v2) n ^ 2 ■ 

61 The open symbols in Figs. E3 and [5] indicate a substantial 

62 difference between the ratios £3/(£2) 3 ^ 2 (a) and £4/(£2) 2 

63 (b) for the MC-Glauber and MC-KLN geometries as in- 

64 dicated. Note as well that the ratios in Fig. [5] are sub- 

65 stantially larger than those in Fig. E3 The latter dif- 

66 ference reflects the different weighting schemes used, i.e. 

67 u>(r±) = r± n and w(rj_) = r± 2 respectively. Interest- 

68 ingly, the ratios for £4/ (£2) 2 imply much larger measured 

69 ratios for ^/(i^) 2 than the value of 0.5 predicted by per- 
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FIG. 6. Same as Fig. [5]for oj(rj 



feet fluid hydrodynamics (without fluctuations) |66|, [67j . 
However, they show qualitative trends which are similar 
to those for the measured ratios tti/faa) < obtained for 
V4 evaluations relative to the ^2 plane [38l [63j. The rel- 
atively steep rise of the ratios in Figs. [5] and [6] (albeit 
steeper for MC-Glaubcr), can be attributed to the larger 
influence that fluctuations have on the higher harmon- 
ics. Note that these are the same fluctuations which give 
rise to the "anomalously low" values of £4 evaluated with 
respect to ^2 in central collisions [f3(| ■ 

Figures [3] - [6] suggests that measurements of the cen- 



trality dependence of the ratios Vs/(v2) 3 ^ 2 and v^/fa) 2 , 
in conjunction with those for V2/V3, W4/W3, v^/v^... may 
provide a robust constraint for the role of initial eccen- 
tricity fluctuations, as well as an additional handle for 
making a distinction between Glauber-like and fKLN- 
like initial geometries. These measurements could also 
lend insight, as well as place important constraints for 
the degree to which a small value of rj/s and/or the 
effects of thermal smearing, modulate the higher order 
flow harmonics [compared to V2] as has been suggested 



23 III. SUMMARY 

24 In summary, we have presented results for the initial 
2s eccentricities £ rl . n <6 for Au+Au collisions with different 

26 weighting schemes, for the two primary models currently 

27 employed for eccentricity estimates at RHIC. The calcu- 

28 lated values of £ n ,n<6; which are expected to influence the 

29 measured flow harmonics v n , suggests that measurements 

30 of the centrality dependence of V2/(«3), "4/^3, «3/(w2) 3 ^ 2 , 

31 V4/ (V2) 2 , etc. could provide stringent constraints for vali- 

32 dating the predicted influence of eccentricity fluctuations 

33 on v n , as well as an important additional handle for mak- 

34 ing a distinction between Glauber-like and fKLN-like ini- 

35 tial geometries. Measurements of v n and their ratios are 

36 now required to exploit these simple tests. 
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